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INTRODUCTION- 

•This  report  presents  a  mathematical  procedure  designed  to  identify 
the  difference  between  two  similar  sets 'of  data  used  in  statistical 
analysis.  It  was  developed  specifically  to  help  solve  detection  prob¬ 
lems  currently  faced  by  the  ASW  (Anti-Submarine  Warfare)  community. 

The  procedure  applies  to  data  whose  statistical  parameters,  such  as 
,  the  means  and  variances,  -aask  essential  information  like  transients 

or  possibly  target  signal  levels  imbedded  in  background  noise. 

Transient  information  in  this  report  implies  that  there  is  a  large  dis¬ 
placement  of  a  very  small  percentage  of  the  data. 

In  this  procedure,  the  variability  in  the  data  is  isolated  and 
defined  by  making  a  comparison  of  the  statistics  of  a  density  function  . 
from  one  time  interval  to  the  next.  This  is  accomplished  primarily  by 
the  manipulation  of  the  differences  between  several  higher  order 
moments.  The  analysis  presented  herein  indicates  that  the  important 
factors  in  this  comparison  are  the  number  of  data  observations  that 
have  been  shifted  in  time  relative  to  a  reference  density  function  and 
the  effective  amplitude.  The  effective  amplitude  will  be  defined  as  a 
displacement  that  will  replace  the  effect  of  the  shift  of  several 
observations  that  are  associated  with  the  transient.  This  type  of  com¬ 
putation  can  be  referred  to  as  a  mathematical  "glimpse"  of  an  approxima¬ 
tion  of  the  variation  in  a  density  function  over  time. 

SUMMARY  OF  RESULTS 

The  report  considers  the  problem  of  extracting  transient  informa¬ 
tion,  that  is,  the  effective  amplitude,  from  large  quantities  of 
ambient  noise  data  without  examining  each  data  sample.  The  effective 
amplitude  is  generally  not  extractable  from  statistical  parameters 
such  as  the  mean  and  the  variance. 

The  task  of  finding  the  effective  amplitude  is  accomplished  by 
j  identifying  the  important  factors  that  describe  the  difference  between 

two  similar  density  functions  derived  from  data  with  overlapping  time 
intervals.  This  difference  is  dependent  on  the  effective  amplitude 
*  and  the  associated  data  samples  within  this  transient.  The  effective 

amplitude  is  then  evaluated  by  the  manipulation  of  the  differences 
found  among  higher  order  moments.  The  knowledge  of  its  value  enables 
the  computation  of  the  number  of  data  samples  associated  with  the 
transient.  Therefore,  it  is  shown  that  knowledge  of  the  factors 
describing  this  difference  will  enable  the  removal  of  transient  informa¬ 
tion  from  large  quantities  of  data  without  the  necessity  of  examining 
each  data  sample. 

CONCLUSIONS 

The  major  conclusions  drawn  from  this  analysis  are  stated  here: 
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1.  This  development  will  enable  the  removal  of  transient 
information  from  large  quantities  of  data  without  the  necessity  of 
examining  each  data  sample. 

2.  This  mathematical  development  is  applicable  to  ocean 
ambient  noise  data  since  these  statistics  display  temporal  fluctuation. 
The  properties  of  the  technique  that  apply  are: 

a.  The  requirement  of  highly  calibrated  measurement 
equipment  is  greatly  re'duced  since  the  analysis  employs  differences 
which  exist  among  means,  variances,  and  higher  order  moments; 

b.  The  information  derived  from  the  analysis  consists  of 
the  effective  amplitude  and  the  associated  number  of  data  samples  that 
change  over  intervals  of  time. 

c.  Application  of  the  mathematical  techniques  presented 
herein  is  directly  applicable  to  the  changing  pattern  of  differences  . 
among  higher  order  moments  of  an  ambient  noise  plus  signal  distribution 
in  a  field  of  sonobuoys  resulting  from  the  effects  of  a  traversing 
target.  Therefore,  this  procedure  affords  a  conceptual  approach  to 
the  detection  of  submarines. 

3.  A  by-product  of  this  analysis  yields  the  interrelations  of 
the  moments  between  two  similar  histograms. 

4.  The  knowledge  of  the  effective  amplitude  and  its  associated 
number  of  data  samples  can  be  used  as  a  measure  for  determining  the 
temporal  fluctuation  between  two  similar  histograms. 
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DISCUSSION 


I .  Background 

Ocean  ambient  noise  is  one  of  the  least  understood  of  the  basic 
parameters  required  to  effective!/  model  acoustic  detection.  One  of 
,  the  difficulties  in  describing  this  phenomenon  is  that  the  statistical 

characteristics  are  from  a  nonstationary  process,  that  is,  the 
statistical  properties  are  time  dependent,  as  observed  in 
references  (a)  and  (b) .  This  effect  can  be  attributed  to  at  least  two 
factors  -  shipping  in  the  ocean  environment,  and  environmental  changes 
in  the  ocean.  Therefore,  the  means  and  variances  of  the  resulting 
distribution  from  these  data  are  not  always  adequate  for  describing 
central  tendencies  and  variabilities. 

Additional  information  contained  within  the  data  can  be  unmasked  by 
considering  the  higher  order  normalized  moments.  For  example,  informa¬ 
tion  about  the  shape  of  the  distribution  can  be  obtained  from  the 
knowledge  of  the  third  and  fourth  moments.  Also,  as  observed  in  refer¬ 
ence  (b),  when  a  transient  condition  exists  in  the  data,  the  higher  order 
moments  become  very  large. 

The  purpose  of  this  report  is  to  introduce  the  development  of  a 
mathematical  technique  that  may  be  applied  to  removing  some  of  the 
difficulties  in  describing  ambient  noise.  -More  important  is  the 
potential  application  of  this  mathematical  development  in  extracting 
transient  information  from  large  quantities  of  ambient  noise  data  without 
examining  each  data  sample. 

II.  Methodology 

This  analysis  is  based  on.  the  premise  that  when  the  comparison  is 
made  between  two  similar  density  functions  with  overlapping  time  inter- 
-.9  vals,  the  resulting  differences  of  the  means,  variances,  and  higher  order 

moments  will  reflect  this  change.  These  differences  are  more  noticeable 
as  the  moments  become  higher.  The  procedure  developed  will  decompose  the 
*  change  into  the  effective  amplitude  and  the  associated  number  of  samples. 

Thus,  the  analysis  determines  whether  the  change  or  shift  in  the  density 
function  is  due  to .a  transient  condition  or  a  small  amplitude  with  a 
large  number  of  samples,  or  some  combination  conditions  between  these 
two  extremes. 

The  basic '.equations  used  in  this  development  will  be  the  sample 
mean,  y 
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X 

Z  Xi 
_  i _ 

N 


i  is  a  running  index  free  1  to  N 
and  N  is  the  niaber  of  data  samples 

Xi  is  the  variate  in  the  data . 

aad  the  equation  for  the  ntn  moment  about  the  mean  is 


v 

n 


N 

l  CXi-u)71 


N 


where  n  is 
if  a  is  equal  to  2,  then  ^2 


the  order  of  the  moment 

2 

is  the  variance  defined  by  the  symbol  a  . 


III.  Analysis  Development  - 
A.  Introduction 


This  analysis  shows  in  detail  how  the  mean  and  variance  are 
altered  by  the  shift  of  one  of  the  variates  of  a  given  data  set.  This 
is  followed  by  a  geometric  interpretation  of  the  shifted  variance  equa¬ 
tion  which  is  then  modified  to  include  the  condition  of  having  more 
than  oae  shifted  data  sample.  The  final  equations  describing  the  shift 
of  the  3rd  and  4tn  moments  will  be  presented,  leaving  the  details  for 
the  appendix.  Manipulation  of  the  shifted  variance  equation  and  of  the 
3rd  aad  4th  monent  equations  will  lead  to  the  solution  for  expressions 
cf  the  effective  amplitude  and  the  associated  number  of  data  points  that 
were  shifted. 

B.  Shifted  Mean 


The  analysis  begins  by  defining  for  an  arbitrary  distribution 
a  reference  mean,  yr 

N 


N 


N 


(1) 


In  a  short  Interval  of  time  one  of  these  observations,.  Xjj  -+.  shifts 
Ajj-j+l  units  while  the  other  samples  maintain  their  original"1  values. 
This  new  observation  is  equal  toofy-j+i  =  %_j+i  +  %-j+l* 
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The  third  and  final  term  becomes  the  difference  between  two  squares, 
since 

2^XK-j»l  ~  Ur  ~  yA^*N-j+l  *  ^H-j+1  _  2^XH-j+l  ~  V^-j+l  *  ^i-j+l 
N  N 

<w  -  va. *  -  qwi  - 

N 


Thus,  the  final  equation  for  the  shifted  variance  is: 

2  =  0  2  +  u  2  +  (XN-j+l  ~  ya  *  ^-jfp2  ~  CXN-j-M  ~  Va)2 

N 


a  =  a  '  +  Va 
a  r  *A 


(7) 


The  variance  for  the  altered  arbitrary  distribution  equals  the  reference 
variance  plus  the  shift  of  the  mean  and  the  term  that  is  the  difference 
between  two  squares.  This  latter  term  will  be  referred  to  as  the  second 

moment  arm.  This  term  will  be  defined  by  the  svmhoi ,  Z«2.  Tt  equals: 


„  2 .  CXN-j+l 
N  = 


~  V, 


,^)2  - 


tXN-j+l 


-»aV 


N 


A  . 

As  N  increases,  decreases  since  an^  tne  second 

i  also  decreases.  Therefore,  as  N  becomes  very  large,  oa2 


moment  arm 

approaches  o^2  which  is  tc  be  expected. 


D.  Implication  of  the  Shifted  Variance 

Before  examining  the  treatment  of  more  than  one  variate  being 
r  shifted,  it  is  worthwhile  to  consider  the  geometric  implication  for  the 

second  moment  arm 

*  „  2  (XN-j+l  ~  ya  +  V-j+l5  "  (XN-j+l  "  va} 

N  =  N 

In  this  discussion,  let  the  data  sample  X^_j+1  be  located  at  the  edge 
of  the  distribution  and  be  shifted  a  distance  Afj-j+l  units  as  indi¬ 
cated  in  figure  1. 
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FIGURE  1.  GEOMETRIC  INTERPRETATION  OF  THE  SHIFTED  VARIANCE 


2 

The  first  term  of  is  equal  to  the  square  of  the  distance 
of  the  shifted  variate  of  Xjj-j+i  from  the  edge  of  the  density  function 
out  to  a  distance  A^-j+i  units,  minus  the  value  of  the  altered  mean  ua; 
that  is,  a  distance  of  (XN_j+i-ua  +  Ajj-j+i)  units.  The  second  term  is  equal  to 
the  square  of  the  distance  from  the  edge  of  the  distribution 
minus  the  distance  of  the  altered  mean  pa;  that  is,  a  distance  of 
(XN-i+l  -  ua)  units.  As  observed,  both  of  these  distances  are  independent 
of  +*e  origin. 


The  moment  arm  associated  with  the  shift  of  the  variate  for 
the  higher  order  moments  becomes: 


n 


(XN-j+l  "  ya  *  AN-j+P  ~  ^N-j+l  ~  ua^ 
N 


The  details  of  the  development  for  this  expression  are  left 
for  appendix  A.  This  particular  term  zft  is  dominant  in  the  computation 
of  the  higher  order  moments .  As  the  order  of  the  moments  increases ,  the 
shift  of  the  distance  AN-j+1  units  from  the  edge  of  the  density  function 

causes  the  moment  arm  to  increase  also. 


Another  term  that  is  important  for  this  analysis  will  be 
referred  to  as  the  second  moment  constant.  It  is  based  on  the  shifted 
variance  equation,  oa^. 
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2  -  a  2  +  v  2  +  (XN-i+1  ~  ya  +  "  ^-j+l  ~  WaJ 


(7) 


N 


or  _ _2-  2 

2  2  2  +  ^N-j+l^N-j+1  11  +  ^N-j+1 

rr  -  rr  -  ll .  =  - - - * - - — - 


°a  -ax  *VA 


N 


Remembering  that  the  shift  of  the  mean  is 


„  =\itl=g 

A  N  a  *r 


The  second  moment  constant,  M-,,  will  defined  to  equal 


2  2  2 

..  °a  -  °r  -  “A  „  .  *»- 

M2  =  - 7 - -  XN-j+l  -  “a  * 

2  y.  J 


2il 


(9) 


Again,  if  X^_j+i  is  at  the  edge  of  a  density  function  and  shifted  out 
to  a  distance  of  A^_^+1  units,  then  the  second  moment  constant  M2  will 
be  the  distance  front  the  altered  mean  out  to  the  distance  halfway 
between  the  edge  of  the  distribution  and  the  transient.  If  the  shift  of 
the  distance  is  inside  the  density  function,  then  Ajj_-i*i  will  be  negative. 
Mathematically,  this  second  moment  constant,  M2,  represents  a  distance 
describing  a  shift  of  a  variate  that  can  be  inside  or  outside  the 
density  function  and  it  equals  the  differences  of  computed  statistics. 

E.  Treatment  of  More  Than  One  Data  Observation 


This  section  discusses  what  happens  to  the  altered  mean 
expression,  ya,  and  the  shifted  variance  equation,  aa  ,  when  more  than 
one  observation  has  been  displaced.  The  genei'al  expression  for  the 
altered  mean  equation  will  be  developed.  However,  the  analysis  considers 
three  observations  being  displaced  in  the  discussion  of  the  shifted 
variance  equation.  It  should  be  understood  that  this  development  can 
apply  in  general  to  many  observations. 

As  shown  by  equation  (3),  the  altered  mean,  ya,  equals 


or 


ua  = 


y  + 
r 


(3) 
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where 


However,  if  there  were  many  observations  being  displaced,  then 


^-i+l  +  ^-1+2 


+  .. .  + 


i=N-i+m 


where  m  is  a  running  index  from  1  to  Ne  and  Ne  is  the 
number  of  data  points  being  displaced 


Let  Ag  be  defined  as  an  amplitude  that  will  replace  the  effect 
of  the  displacement  of  many  observations.  It  equals 


Ag  =  jj—  I  Ai 


Now,  ua  will  equal 


e  i=N-j+m 


N  A 
e  e 


and  Ag  will  be  referred  to  as  an  effective  amplitude.  It  will  be  deter¬ 
mined  by  the  use  of  the  differences  of  the  higher  order  moments.  The 
value  of  Ne  is  determined  by  the  knowledge  of  ua  which  is  equal  to  the 
difference  between  the  altered  mean  and  the  referenced  mean.  Therefore,  Ne 
is  defined  as  the  effective  number  of  data  points  that  are  associated 
with  the  shift  of  the  mean  va>  ° r  the  number  of  data  points  that 
significantly  contribute  to  the  difference  between  moments  of  similar 
statistics. 

The  modification  of  the  shifted  variance  equation  for  the 
displacement  of  three  observations  will  be  developed  by  first  referring 
to  equation  (5).  It  should  be  observed  that  the  reference  variance, 

°x2,  and  the  shift  of  the  mean,  ua>  will  not  be  modified  when  considering  the 
shift  of  more  than  one  data  point.  The  moment  arm  expression  is  the 
term  that  will  have  to  be  modified  and  the  form  used  (page  5,  second  term  of 
equation  (6))  in  this  discussion  is: 
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2AN-H 


»N-1 


.  ,  -  u  + 
3+1  a 


As  discussed  earlier,  the  expression  in  the  bracket,  XN  -  y  pa  +-  -4-* 1 
geometrically  represents  a  measured  distance  halfway  between  the  variate 
before  being  shifted  and  to  the  distance  of  the  transient.  If  there  were 
several  observations  shifted  approximately  A  units,  then  with  each 
observation,  there  would  be  an  associated  distance  shift  of  the  second 
moment  arm.  Since  considerations  are  now  given  to  data  with  more 
than  one  variate  being  shifted,  it  is  worthwhile  to  consider  another 
term  to  be  defined  as  Xe.,  the  effective  distance.  It  is  equal  to 
the  average  value  before  the  shifting  of  all  variates  that  are  shifted. 
The  techniques  to  be  developed  are  the  approximate  value  of  the  effec¬ 
tive  amplitude,  Ae,  that  is  related  to  the  shift,  and  the  associated 
number  of  data  points,  Ne,  associated  with  the  shift. 

The  modification  of  the  second  moment  arm,  Zj^2,  due  to  the 
shift  of  few  variates  that  were  caused  by  amplitudes ,  Aj.  .  . ,  Aj.  .  _, 
a  is  shown  below:  ~^+ 


331(1  Vj+3 


2  2  - 
ZN 


2AN-;h 


~  ^a  + 


Vi+J 


2AN-j+2  fv 

N  XN-j+  2~ 


^XN-j+  3  "  Ua  + 


The  terms  in  the  brackets  of  the  expression  for  are  also 
associated  with  distances  that  are  shown  below  in  figure  2. 


Xe  -  +  f 


,  .liVi*1  ^ 


Vj+3 

AN-3'+2 

^  \l 


nTti 
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The  technique  developed  in  this  section  shows  how  to  solve  for 
an  average  or  an  effective  amplitude,  Ae.  Once  Ae  is  known,  Ne  can  be 
computed  by  the  knowledge  ofy^.  The  additional  term,  Xe,  can  also  be 
evaluated.  This  distance  Xe  will  normally  be  within  the  density  function. 
Figure  2  also  shows  the  relative  location  of  Xe  as  well  as  Ae/2. 

Since  this  computation  considers  an  effective  amplitude,  Ae, 
and  an  effective  distance,  X£,  errors  resulting  from  the  approximation 
exist  in  the  analysis.  However,  in  the  case  of  a  transient,  where  the 
shifts  are  large,  the  approximation  error  should  be  very  small. 


In  order  to  solve  for  the  terms  Xe  and  Ae,  they  have  to  be 
introduced  within  the  brackets  of  equation  (12) .  This  can  be  accom¬ 
plished  by  setting  each  bracket  term  plus  a  differential  equal  to 
(Xe-Ua+Ae/2)  as  shown  below. 


fXe  '  ya  +  2  1  "  ^-j+l  "  ya  +  ^  +  6N-j 


5+1 


(13) 


”  ^XN-i+2  '  ua  +  ~2_1  +  6N-j 


-5+2 


5+2 


-XN-j+3  “  pa  +  2  ^  +  6N-j+3 


The  differential  term  will  adjust  the  X's  and  the  A's  to  be  equal  to  Xe 
and  Ae,  respectively.  Thus,  introducing  an  error  into  the  bracket  terms, 
and  rewriting  equation  (12),  we  obtain: 


7  2  2AN-3+l  rY  S  k  1 

ZN  "  N  fXe“  Pa  2  “  6N-j+r 


N 


[X£ 


^a  +  f 


‘  6N-j+2^ 


(14) 


+ 


2A, 


N-j+3 

N 


6N-j+3 


] 


2  ^AN- j+1  *  Vj+2  *  Vj+3^ 
N 


N  (AN-j  +  l  6N-j+l  +  ^+2  6N-j+2  +  V-j+S  5N-j+3^ 
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Since 


VA  = 


K 

e  e 


N 


equation  (12)  becomes: 


^-j+l  +  ^N-j+2  *  ^N-j+3 
N 


2N  A 
e  e 

N 


t(Xe-  ».♦£>-  'Wf 1  *  ‘"tr  ^  05) 

*N-j+l  +  Vj+2  Vj*3 


where  the  S's  can  be  both  positive  and  negative  values.  The  errors  for 
the  second  moment  arm  should  be  very  small  and  therefore  are  ignored. 


The  second  moment  constant  in  the  case  of  several  shifts  would 
be  equal  to: 


Thus,  there  is  one  equation  with  two  unknowns. 

In  order  to  solve  for  Xe  and  Ae,  additional  equations  are 
needed.  They  are  found  by  introducing  a  shift  of  a  variate  into  the 
higher  order  moment  equations.  The 'third  and  fourth  higher  order  shifted 
moment  equations  are  developed  in  appendix  B  and  are  given  below. 

The  third  moment  shift  equation: 

M3a  ’  u3r  -  3Vr  ‘  "a3  *  ZN 3 

The  fourth  moment  shift  equation: 

*J4a  =  y4r  ‘  4yA  w3r  +  6yA2  ar2  +  V  +  ZN4  ™ 

and 

Z/  =  r  [(Xe  "  ua  +  Ae)n  -  (Xe  -  Va)n] 

As  mentioned  above,  the  approximation  errors  in  the  above  equations  are 
considered  negligible.  The  subscripts  (a)  and  (r)  refer  to  the  altered 
distribution  and  the  reference  distribution,  respectively. 

The  third  and  fourth  moment  constants,  which  are  developed  in 
appendix  B,  are  given  below.  The  third  moment  constant,  neglecting  the 
approximation  error. 
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„  p3a  p3r  +  3yA  CTr  ‘  yA  rv  Ae^?  .  1  ,.Ae^ 

M3  =  - — - -  CXe  -  Ua  +  —)2  +  3  C— )2 

3ua 


(19) 


For  the  fourth  moment  constant. 


M4  = 


V_U  -  |14r  *  4uA  "Sr  ~  6mA2  °r2  ~  V 

*  ' 


(20) 


Ae  Ae  Ag 

=  (Xe  -  Pa  *  y-)  [(Xe  -  ,ia  +  2~)^  +  02) 

Thus,  two  more  moment  constants  relating  Xe  and  Ae  exist. 

It  is  observed  that  since  two  unknowns  exist,  there  are  more  equations 
than  needed.  Indications  are  that  the  higher  the  order  cf  the  moment 
constants  used,  the  more  accurate  the  answer  for  Xe  and  Ae. 

F.  Solutions  for  the  Differences 


Solving  for  Xe,  Ae  and  Ne  using  the  constants  M2  and  M3,  the 
second  moment  constant  becomes 


a  2  _  0  2  _  y  2 
a  r  A  „  e 

H,  - -  =  Xe  -  ua  + 

2uA  • 

For  the  third  moment  constant: 


'a  '  2 


+  3wA  °r  +  yA  Ae  2  1  Ae  2 

M3  =  ^ - —  =  (Xe  -  ya  +  ^-)  +  i  (/) 

3p. 


m3=  oy2 

Solving  for  Ae 

Ae  =  '12  (M3  -  M22) 

Hence , 

Ne  =  —A.  and  Xe  =  M2  +  Fa  -  Y~ 

where  Ne  is  the  associated  number  of  samples  and  Xe  =  the  relative  position 
of  the  variate  before  the  amplitude  of  Ae  units  is  added  to  the  data. 
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Solving  for  Ae  using  M2,  M3  and  M4 

M4  =  y4a  ~  P4r  *  %  u3r~  6pA2  gr2  '  V 

4p* 


ft*.  -  -  r’2  *  tr-)2i 

However,  the  bracket  term  is: 

[O'®  -  “a +  r-)2  -  'r-)2]  -  [«e  -  »a  *  r-)2  -  i  4%2  *  1 1^.)2) 

2  A- 

-  m3  *  I  Or’2 

The  fourth  moment  constant  equals 


M4  -  (Xe  -  ua  ♦  ^)  [M3  *  1  4%*] 


*  «2  [Hj  ♦  |  fcpl2) 


Solving  again  for  Ae, 


Ae  =  ^6 (M4/M2  -  M3) 

It  is  observed  that  M4,  M3,  and  M2  are  functions  of  the 
difference  of  the  computed  moments  and  are  not  dependent  on  absolute 
measured  values.  Both  Ae  and  Ne  are  dependent  on  the  moment  constants 
ana  independent  of  absolute  values.  They  are  the  common  factors  that 
interrelate  two  similar  sets  of  statistics.  Therefore,  they  identify 
the  difference  between  two  overlapping  density  functions.  These 
ingredients  that  describe  the  differences  can  be  an  aid  in  the  analysis 
cx  time  varying  statistical  data  such  as  ocean  ambient  noise. 
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IV.  Gomdcsioos  of  the  .Analysis 


Techniques  developed  from  this  analysis  yield  a  new  tool  that 
extracts  ■fee  differences  between  two  sets  of  similar  statistics.  The 
extraction  will  be  the  displacement  (A-).  that  will  replace  the  effect 
of  a  shift  of  Several  observations,  and  fee  number  of '  observations 
(S*).  These  two  quantities  reflect  fee  changes  feat  alter  statistics 
feat  are  measured  from  one  tine  interval  to  the  next.  The  analysis 
models  fee  occurrence  of  these  effects  even  if  they  are  deeply 
imbedded  in  fee  temporal  ocean  arbient  noise  -data.  The  important 
feature  of  this  tool  is  that  the  gross  changes  in  the  data  can  be 
resDorved  without  the  necessity  of  examining  every  data  sample. 

Consider  the  situation  of  ocean  ambient  noise  levels  being 
monitored  by  the  field  of  soncbuoys.  If  a  target  traverses  this  field, 
fee  tine  varying  averages  of  fee  ambient  noise  levels  may  mask  informa¬ 
tion  relating  to  its  presence.  However,  the  changing  means,  variances, 
and  sene  of  the  differences  of  the  higher  order  moments  will  reflect 
fee  aunoer  of  data  samples,  Ne,  feat  vary  in  time  throughout  the  field. 
Possible  target  presence  can  be  indicated  by  an  increasing  value  of  Ne 
relative  to  fee  data  of  the  other  buoys.  If  a  transient  were  present, 
fee  value  of  Ag  would  increase.  These  Measurements  are  dependent  on  the 
relative  changes  of  moments  and  not  the  absolute  value  of  highly  cali¬ 
brated  buoys  at  each  position  in  the  field.  Therefore,  the  use  of  this 
mathematical  technique  affords  a  conceptual  approach  for  the  detection 
of  submarine  targets - 

The  tool  can  be  used  in  aiding  fee  analysis  of  temporal  fluctuations 
associated  with  ocean  ambient  noise  data.  A  by-product  of  this  analysis 
yields  interrelations  of  moments  between  similar  density  functions. 
Therefore,  fee  tool  can  aid  in  measuring  fee  departures  of  time  varying 
statistics. 
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APPENDIX  A 

THE  SHIFTED  HIGHER  ORDER  MOMENTS 


In  the  development  of  equation  C7) »  the  shifted  variance  equation 
,  is  expressed  as 


a  2  =  a  2  +  p.2  + 
a  r  A 


y 


(7) 


It  is’  observed  that  oa2  is  a  function  of  the  reference  variance,  a  term 
due  to  the  shift  of  the  mean  and  the  moment  arm.  It  will  also  be  shown 
that  the  higher  order  moment  equations  will  have  the  same  form. 

In  order  to  simplify  the  development  of  the  higher  order  moment 
equations,  the  Kronecker  delta  will  be  introduced  into  this  section  of 
the  paper.  This  term  has  the  following  properties. 


5.  .  =/1  1  = 

l,N'3  \o  i  i  N-j 


Now,  let  the  variate  Xi  be  equal  to 

Xi  =  Xi  + 

Hence,  when  the  summation  of  Xi  from  i  to  N  is  performed,  the  amplitude  of 
A  units  that  causes  the  shift  in  variate  will  only  have  a  value  when 
i  =  N- j . 

Development  of  the  Shifted  Third  Moment  Equation 

Remembering  that  Pa  =  Pr  +  PA  anc*  starting  with  the  third  moment 
equation 

1  N  i  N 

y3a  =  N  E  (Xi"pa)3  =  jf  Z  (*i+  6i,N-j  AN-j  -  Pr^3  (A-1) 
i  i 

Let  Zi  =  Xj-pr-pA  and  =  Xi~pr 

Now  consider  the  third  moment  expression  when  i=N-j,  then 

y3aji=N_.  =  S’  (zN-j  *  3 

“  FT  H-j  +  3AN-j  ZN-j  +  4-j  ZN-/+ 

=  w  +  ^-j^N-j-V^  +  +  Vj] 


A-l 
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The  tens  associated  with  the  moment  arm  will  be 


y  ■  k  '3AN.j « N-j-y 2  *  3an-j  ®shV  *  i 

Adding  and  subtracting  yields  (X^  .  -  3 

y  =  N  *  3AN-j®N-j'11a'12  +  3AN- j (XN-j "'“a5  (A‘3) 

*  *N-j  -  V)  '  "a)3] 

ty-i^a  -  y3 

N 

Thus,  the  third  moment  expression  when  i  =  N-j  is 


u3a  .  “  M  <Xi-V‘,A)  *  2 

l=N-j 


3  j.  7  3 
N 


*  S’  Cyi'-“A)3  *  ZN3 

Now,  consider  the  expression  cf  u^a  with  the  amplitude  A^_-  term  being 
removed  from  the  summation 

"3a  =  S’  ?  ^‘-"A)3  *  ZN3 


N  3p.  N 

1  —  3  _  _A  5  YZ  + 


3"a  ; ..  %3 . ,  3 


_  ±  y  v  3 _ "  r  v  z  +  _ —  T.  v _ : 

~  N  .  yi  N  .*  yi  N  .  yi  N 
i  x  .1 


+  Z, 


N 


■  "3r  '  3"a  °r  -  "a3  *  %3 


The  summation  of 


N  N 

N  ?  yi  =  N  l  CXi”pP  =  0 


A-2 
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\ 


Therefore,  the  third  moment  equation  due  to  a  shift  in  the  variate  by 
A  units  is 


"3a  =  "3r  '  3"a  “r*  '  "A3  *  V  tA-4) 

Considering  the  case  when  more  than  one  variate  has  been  shifted  but 
ignoring  the  approximation  errors,  becomes 

=  N~  +  ^e)3  ~  $e  ~  Pa)3]  (A-5) 


Development  of  the  Shifted  Fourth  Moment  Equation 
Starting  with  the  moment  equation 


"4a  =  W  3  CCi-"a5'*  *  k  3  0[i-"r-l'A+«IH  V.)" 
but  first  by  just  considering  the  term,  p4&,  when  i  =  N-j,  we  obtain 


N 


4a 


i=N-j  *  N  i  tZN-j  *  M5 


N  tZH-j  *  4AN-j  zN-j  *  6AN-j  zN-j  *  4AN-j  zN-j  *  flN-j) 
The  moment  arm  part  of  the  expression  is 

V  =  N  (4AN-j  zN-j  +  Vj  zN-j  +  4AN-j  zN-j  + 

=  i  (zjj.j  +  4A  z*_.  +  oA^  Z2_.  +  4A~  zN_.  +  A4  -  z^_.) 
=  (ZN-j  *  "  *N-j 


N 


7  4  _  ^N-j  wa  +  AN-j}  (XN-j  "  V  •  z 
4n  -  - t - i - l - »  zm^  - 

N 


N-j  “  ^-j  -  ya  (A-6) 


A-3 
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Considering  the  case  when  more  than  one  variate  has  been  shifted  but 

ignoring  the  approximation  errors  for  Z  4 

N 


N 

V  =  ir  f(Xe  -  Va  +  V4  -  exe  “  Va)*1]  (A-7) 


Thus,  the  fourth  moment  expression  when  i  -  N-j  becomes 


4.i 


i-N-j 


-w«i- 


V 


Z  4 


V>4  +  *N 


4 


j 


«• 


Now,  consider  the  expansion  of  u^a  with  the  amplitude  Ajj  -  not  being 
considered  in  the  summation 

O'!'  V*  V 

1 


N 


1"  u  %!!  i  V?  ,  4pA3^  i  NyA 

=  N  x  ^i  -  IT  1  yi  +  IT-  1  yi  -  —  £  ^i  +  ~ 


v4a  =  p4r  '  4vk  v3r  +  S*  °f  +  "a  +  V 


(A- 8) 


The  term 


N 


Z  yi  =  N  Z  (Xi_lJr)  =  0 


As  can  be  seen  in  this  development,  the  expression  for  the  n^  order 
moment  of  the  shifted  variate  is 


“na  ■  if  E  (yi-»A)n 
1 


+  ZNn;  where  =  Xi-yr 


(A-9) 
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APPENDIX  B 
THE  MOMENT  CONSTANTS 


The  third  moment  constant  will  be  developed  by  rewriting 
equation  (A-4)  and  ignoring  the  approximation 

y3a  "  y3r  +  3yA  °r2  +  PA*  '  Ti  ^Xe  '  ya  +  Ae)3_  (Xe  ~  ya)3l  C8'1) 

N 

y3a  '  y3r  +  3yA  +  yA 3  =  iT  t3AeCxe  -  Pa)2  +  3Ae2(xe  -  Pa)  +  Ae3] 


3N  A  A  2 

-g-5-  UXe  -  Pa)2  *  AeCXe  -  Pa)  *  -f-J 


3N  A  A  2  A  2  A  2 

e  e  «• e->  e  +  e  1 

4  3  1 


=  CCXe  -  Pa  +  r> 


3N  A 


P*„.  -  P*„  + 


2P/v  0  2  +  P43  =  -  %  6  UXe  -  Pa  +  +  T 


3a  3r  A  r 


N 


a  2 


3  l2 


The  third  moment  constant  becomes 


..  “3a  '  w3r  +  3llA  *  vk  .  ..  .  V,  .  1  \2 

M„  —  1  -r  »i  «  fat  —  (Xg  —  Pa  "*■  9  )  T  Co  ) 


3  N  A  /N 
e  e' 


But 


M  A 

e  e  , 

"a'TS—  and 


,,  11 3a  ~  +  3uA  °r  UA  fv  Ae,  2  .  1  2 

H3  *  - 3T; - '  Xe  -  “a  *  r 5  *  *  ’ 


3  '■  •2 


The  fourth  moment  constant  will  be  developed  by  rewriting 
equation  (A- 7)  and  ignoring  the  approximation 


B-l 
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V  =  r-  t«e  -  »a  *  *•)"  -  CXe  -  i.a)*i 


N, 


jp  [4Ae(Xe  -  ya)3  +  6Aq2  (Xe  -  ya)2+  4Ae3  (X6  -  Ug)  +  Ae4] 


4NeA6 


=  HP  [(Xe  -  ya)3  +  f  Ae(Xe  -  ya)2  +  Ae2  (Xe  -  ya)  *  4-  1 


4NeAe 


jp  [(Xe  -  ya)3  +  3(^)  (Xe  -  ya)2  +  3(^.)2  (Xe  -  ya) 

3 


i  3  .2 

1  Ae  Ag 


♦f-f— 


4N  A  A  3  A.  2  A 

HP  [(Xe  -  Pa  ♦  p)  ♦  (p)  (Xe  -  ya  *  ~)] 


4N  A  A  A  .  A  „ 

V  -  -F  Cxe  -  “a  +  F>  CCXe  -  »a  *  Y>  +  0  1 


(B— 3) 


Rewriting  equation  (A-8)  and  substituting  equation  (B-3)  into  it  yields 


“4a  '  “4r  *  4“a  “3r  '  S*  °r  '  “a"  * 


4N  A 
e.  e 


iP  <Xe  '  »a  +  t(Xe  "  ua  +  +  1 


Dividing  both  sides  of  the  equation  by  4NeAe/N,  the  fourth  moment  constant  is 
obtained. 


M  ~  ‘‘4a  '  U4A  I  4"a  t‘3r  '  6llX2  °r2  ~  “a" 

4  4  NeAe/N 

A_  A  2  A  2 

“  <Xe  -  *a  +  7*  f(Xe  "  ^a  +  7*  +  3 


(B-4) 


B-2 


